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We present a novel and simple mechanism to produce vortices at any desired spatial locations in 
harmonically trapped Bose-Einstein condensates (EEC) with multicomponent spin states coupled 
to external transverse magnetic fields. The vortices appear at the spatial points where the spin- field 
interaction vanishes and, depending on the multipolar magnetic field order, the vortices acquire a 
corresponding charge. We explicitly demonstrate our findings, both numerically and analitically, 
by analyzing 2- and 3D BEC via the Gross-Pitaevskii equation for atomic systems with either two 
or three internal states. We further show that, by an spontaneous symmetry breaking mechanism, 
vortices can appear in any spin component, unless symmetry is externally broken at the outset. 
We propose an experimental realization of this scenario using an ultracold gas of ^^Rb occupying 
all three F — 1 states in an optical trap. We suggest that a further time varying of the external 
magnetic field may lead to a route for quantum turbulence. 

PACS numbers: 67.85.Fg, 03.75.Lm,05.30. Jp 



I. INTRODUCTION 

Quantum vortices are the hallmark of superfuidity in 
either fermionic or bosonic ultracold gases [IjiT^. Follow- 
ing the seminal work of Abrikosov^ and other works on 
vortices in rotating Helium [9l [10], it has been widely 
verified that rotating Bose-Einstein condensates (BEC), 
under appropriate angular velocities of rotation [SHTl [TT] , 
indeed show a triangular vortex lattice as the most stable 
state. Alternative ways have been devised to create vor- 
tices, such as phase imprinting [T2]-[15] and, very recently, 
by creation of artificial gauge fields [TG^.TT. In this article 
we introduce another way of creating vortices; one that 
necessarily requires the presence of more than one atomic 
internal spin state coupled to an external transverse mag- 
netic field. In this regard it is of relevance to recall 
that, even though very early in the development of opti- 
cal traps [m [19] it was pointed out [20 -22 the additional 
richness of spinor or multicomponent Bose-Einstein con- 
densates, most of the efforts in creating and studying vor- 
tices were focused on rotating condensates [23 . The nov- 
elty of our proposal, in addition to exploit the spinorial 
structure of the condensates, is that the vortices can be 
created essentially at any desired spatial location. This 
property, as we shall see, depends on the ability to pro- 
duce a corresponding on-demand tailored magnetic field. 

The essence of vortex appearance follows from the 
same argument provided by Abrikosov[8 . In supercon- 
ductors, vortices appear in the spatial points where the 
type II superconductor allows the magnetic field to pen- 
etrate. In those vortices (or flux tubes as they are also 
known) the density is zero and the superconductor cur- 
rent circulates around the vortex to cancel the field in the 
bulk. In our proposal the situation is very similar, except 
that we work with multi-component wavefunctions and. 
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unlike Abrikosov, we work with non-uniform magnetic 
fields. 

We find, interestingly, that the vortices are nucleated 
not where the magnetic field is able to penetrate the con- 
densate, but on the spatial points where the field is zero. 
Additionally, the vortex solution will not have the same 
charge in all spin components of the wave function. That 
is, if a vortex of charge Q is nucleated on one spin com- 
ponent, the vortex solution on the other spin components 
for the same spatial location must have a charge q Q. 
Moreover, if this charge is zero, ^ = 0, no vortex will ap- 
pear and the spin component will develop a density accu- 
mulation or spike so as to cancel the vortex formed on the 
other component (s), leaving the total density unaltered 
in the usual gaussian-like form. In this way, vortices will 
be nucleated "on demand" at the spatial positions where 
the magnetic field is zero, albeit not with the same charge 
on all spin components. 

This approach presents some additional interesting 
consequences. The first one is that the vortices can be 
created with arbitrary integer charge, Q = 1,2,3..., by 
selecting a field with appropriate multipolar order; dipo- 
lar fields give rise to Q = 1, quadrupolar ones to Q = 2, 
et cetera. In this way, not only the position of the vor- 
tices, but their charge as well, can be controlled using 
on-demand tailored magnetic fields. 

The second interesting consequence of our approach 
is that it is possible to assemble vortex lattices "on de- 
mand" by selecting magnetic fields with zeroes on the 
lattice sites. Now, since all solutions obtained from the 
Gross-Pitaevskii and Ginzburg-Landau equations are of 
minimum free energy, we conclude that triangular vortex 
lattices observed both at superconductors [24 and BEC 
alike [5]-[7l |TT] , are a peculiarity consequence of the use of 
uniform magnetic fields. Rotating BEC are not an excep- 
tion to this. It is possible to see that a rotating BEC is 
equivalent to a superfluid in a uniform magnetic field. In 
this way, our approach may have the further advantage 
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of providing evidence of a unique cause underlying the 
different vortex nucleation techniques currently at use. 

In this article we first analyze the corresponding Gross- 
Pitaevskii (GP) equation to describe how the vortices ap- 
pear. We consider a very general case and then illustrate 
several different cases in 2 and 3D gases by numerically 
solving the corresponding GP equation. In the last sec- 
tion we suggest how this situation can be experimentally 
realized in a gas of ^^Rb occupying all three F = 1 states 
in an optical trap. In this regard, we point out that pre- 
cursors of the present study are the analysis of Refs. [25] 
and [26] , the former discussing the appearance of vortices 
in the center of a loffe-Pritchard trap in a F = 1 conden- 
sate, while the latter being a detailed study of vortices 
on two- level condensates, yet, in rotating systems. 



II. A MULTICOMPONENT BEC IN A 

TRANSVERSAL MAGNETIC FIELD 

Consider a 2 or 3 dimensional BEC superfluid at zero 
temperature. The gas is confined by an optical trap. 
In addition, there exists an external transveral magnetic 
field that couples to the x and y spin components. The 
atoms have spin J. The 2 J + 1 GP equations are. 
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where r is the D-dimensional position vector. The Greek 
letters, a,/3 = l,2,...,2J+l, denote spin components. 
The confining external potentials are V^^{f) and we have 
considered the possibility that each spin component feels 
a different confining trapping potential; however, for il- 
lustration purposes we shall assume below that they are 
all the same harmonic potential. In general, we assume 
the external trapping potentials to be smooth at any 
point. We can also, in principle, assume that the scat- 
tering length is different for any pair of components; but 
again, in the numerical evaluations below, we will con- 
sider the simple case = This is certainly an 
important point since it is known[27 that the scatter- 
ing length truly depends on the internal atomic states 
and this in turn leads to different types of macroscopic 
"magnetic" properties of the condensates. By assuming 
a unique positive scattering length throughout, one may 
assert that the condensates with which we are dealing are 
at the border of being ferromagnetic or polar [20 . This 
does not invalidate the point we are highlighting here. 

The most important part, and the novelty of our argu- 
ment, is the coupling of the spin components via an ex- 
ternal inhomogeneous transversal magnetic field B{f) = 
Bx{T)ex + By{f)ey. Accordingly, J^^ and Jy^ are the x 
and y angular momentum matrices of spin J. In addition, 
we also include a uniform longitudinal magnetic field Bz . 



We consider here the cases J = 1/2 and J = 1 only, but 
the argument may be extended to higher orders. In ([2|, 
/i is the chemical potential and mo the magnetic moment 
of the atom. 

Now, the external magnetic field must obey V • ^ = 
at any point. In addition, since we want to resemble 
a true experimental situation, we should demand that 
throughout the condensate V x 5 = is also obeyed. 
In what follows we will consider these two restrictions to 
build the proposed fields, however, we will also include 
more general situations. We also defer to the following 
sections the question of how the proposed fields may be 
produced. 

The main point of the discussion is based on proposing 
that the transversal magnetic field becomes zero at iso- 
lated points fo, namely, Bx{fQ) = By{fQ) = 0. Placing 
the origin at fo, the magnetic field very near that point 
is assumed to have the following form, 

B^x, y) « (l^-'(7r/2, cj>) + (-l)'l^'(7r/2, cj>)) (2) 

and 

By{x, y) ^ -iBi {Yf\n/2, 4>) - {-l)%\n/2, 4>)) r' (3) 

where are the usual Spherical Harmonics of order 

/ = 1, 2, 3 . . ., tan^ = y/x and = -f 7/^. Because the 
fields are constant in the z direction, the fields are to be 
evaluated at = 7r/2. We shall call / = 1 dipolar, / = 2 
quadrupolar, etc. [28]. Bi are constants with appropriate 
units. We do not write the constant z component of the 
field. For the sake of exemplifying the two simplest cases, 
we write down the explicit forms for / = 1 and / = 2, 



B{x, y) ^ Bi {xe^ - ycy) for / = 1 



(4) 



and 



B{x,y)^B2{{x^ -y^)ex-2xyey) for 1 = 2. (5) 

Let us turn now to the main problem: vortex nucle- 
ation using the fields already discussed. For simplicity 
we analyze the case of spin J = 1/2. This case may be 
realized with the coupling of the magnetic field to only 
two components of a larger actual bosonic state, F > J 
with F integer. In this case, Jx and Jy are the x and y 
Pauli matrices. 

Let us analyze equation ([2| in the vicinity of fo. In 
a spinor notation, we denote a = 1 as and a = 
2 as ". Hence, Q reduces to the following coupled 
differential equations. 



V' + KUr^ ro) + g {\^^{r)f + \^_{r)f) 



2m 

^moBir'e^''^^{r) = /i±^±(r) 



(6) 



where /i± = /i ± mo^^- Since the external potential is 
smooth, we can set Vg^^(r + fo) ~ Kttl^o)^ namely a 
constant near fo. In such a vicinity, we can fairly see 



3 



that there exists a solution with cyhndrical symmetry of 
the form, 



as, 



^+(r) « /|(r)e^'^ and « /i(r) 



(7) 



or, alternatively, with the solution without a phase 
and with the " with the phase exp(— 

Substituting the above forms into Eqs.(|6| yields two 
coupled equations for the functions /+(r) and /i(r). The 
ensuing equations can have a solution only if near r = 
they behave as, 

/|(r) ^ C+r^ and f_{r) ^ as r ^ 0. (8) 

This can be seen by proposing a formal power expansion 
near r = of both functions [8]. 

The component '0+(r) represents a vortex solution of 
order (or charge) /. This is readily seen by calculating its 
circulation 



Vo - dr 



n 



m 



-27r/, 



(9) 



where Vs is the superfluid velocity field and the integra- 
tion is on any closed contour around tq. 

The component ip-(f) is not a vortex solution; namely, 
its circulation vanishes C_ = 0. However, it can be shown 
that its value at the vortex location fo obeys. 



1 



\^-{ro)\' = — {^i-V-,{T,)) 



(10) 



which is the Thomas- Fermi solution at that point. This 
is the density accumulation or spike referred to at the 
introduction. 

Had we chosen the solution with the phase in the " 
component, a vortex would have been obtained in that 
component with negative circulation, C_ = —{h/m)27rl. 
Analogously, C+ = and no vortex would have been 
formed for the component. 

Within the previous discussion, it is important to note 
that there is no a priori argument to decide where the 
vortex will be formed. A BEC with the same potentials, 
= V~ can nucleate vortices in either component de- 
pending on the initial conditions. There must be a spon- 
taneous symmetry breaking mechanism making the de- 
cision. As we shall see below, the constant longitudinal 
field Bz breaks such a symmetry and the vortices do ap- 
pear in the positive direction of the field. 

Now, let us note that since the above argument only 
gives the solution at the given point fo where the external 
magnetic field vanishes, one can consider cases where a 
magnetic field has zeroes at different spatial locations. 
Applying the same argument near any of those points, 
a vortex will be nucleated at any desired position. It is 
also a simple calculation to verify that the "size" ^ of the 
vortex is the same one as that of the classical solution of 
Gross[29 and Pitaesvkii[30l, namely ^ {Ti^ /2miiY^'^ . 

The extension to J = 1 is immediate. To see this, 
consider the magnetic field term only. It can be written 
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This yields a vortex at Mp = +1, no vortex at Mp = 
and a — / vortex at Mp = — 1. The opposite solution is 
also equally valid. 

At this point it is relevant to comment on other vortex- 
like solutions to the GP equation. Because of the non- 
linearity of the GP equation, other solutions besides the 
one obtained by variable separation may be valid. Until 
now we have considered only one kind of solution in both 
the J = 1/2 and J = 1 cases. In the former, all the 
topological charge is concentrated in only one of the spin 
components or "-" ; in the later, equal and opposite 
charge is found in the and spin components. 

However, these are not the only valid solutions to the 
GP equation presently considered. Let us examine the 
J = 1/2 and J = 1 cases. 

In the J = 1/2 case, close to the vortex where the 
solution behaves as in eq.([7|), it can be seen by inspection 
that the GP equation accepts solutions of the form 



iP+{f) ^ /|(r)e'°* and ^/;_(r) ^ /i(r)e 



lb4, 



(12) 



where a and b are integers and a = b -\- 1. In this way, 
there are / + 1 solutions to the GP equation presently 
considered. 

In the J = 1 case, the rule to select possible solutions 
is slightly different. Using similar arguments to those 
leading to eq.(12), it can be shown that solutions of the 
form 

^+(0 « fi{r)e'-f,Mn « /^(r)e-*,^_(r-0 « f_{r)e''"f 

(13) 

where a,b and c are integers, a = c-\-l and b = c — I, are 
also valid solutions. In this way, in the J = 1 case under 
an multipole field, there are / + 1 solutions as well. 



III. A GALLERY OF VORTICES ON DEMAND 

In this section we illustrate the present mechanism for 
vortex formation with a variety of cases in both 2D and 
3D; additional details case by case are also provided. For 
simplicity, several assumptions are made. The trapping 
potential is assumed to be an isotropic harmonic poten- 
tial for all spin components with the same frequency cu. 
In all cases the same scattering length is assumed for all 
collisions, namely, = gSa(3^ with g a free parameter. 
Units in which h = m = uo = 1 are used and hzi = m^Bi 
is defined as the parameter that determines the strength 
of the external magnetic field. Throughout, the field in 
the z component is assumed constant, Bz = BiZq. In all 



cases below, the corresponding set of GP-equations is nu- 
merically solved with the method proposed by Zeng and 
Zhang [3T]. The chemical potential is part of the solution 
and it is also our convergence parameter; its value with 
its corresponding error is given in all cases shown. The 
density plots as well as phase- velocity fields Vs = V<l>, 
with ^ the corresponding phase of the tpa{^) component, 
are also shown; the circulation is plotted as well. 



A. Single / = 1 and / = 2 vortices in 2 and 3D BEC 

This is the simplest case to consider and we begin with 
two internal states, namely J = 1/2. Letting / = 1 in 
Eqs.(2| and ([3|, the field is dipolar in the sense that it 
is a linear combination of spherical harmonics of order 
one. It is very similar to the loffe-Pritchard trap field 
except that our field has a constant z component and, 
consequently, the sign of the y component must change 
to preserve Maxwell's equations. 

Fig. [l] shows a vortex at the origin fo = in a 2D 
BEC and Fig. [2] a vortex at an arbitrary fo, see figure 
caption for details. As mentioned above, the density and 
phase-velocity plots, as well as numerical calculation of 
circulations, are also shown. Square contours around the 
vortex were used for numerical integration of the circu- 
lation; as the contour includes enough points, one of the 
circulations gives C/27r = ±1 while the other gives zero 
within numerical accuracy. 

The population of both states, p± = |'0±(r)P and the 
total one, p+ + p_, are also indicated in the figures. The 
values for both the and " states are close to 0.5 
but not exactly. We have no explanation for these devia- 
tions from half-and-half population. It can also be noted 
from the figures that if a vortex is nucleated on one of the 
spin components, no vortex will be nucleated on the same 
spatial position in the other components. Additionally, 
the other spin components develop a density accumula- 
tion or spike that cancel the vortex formed in the other 
component, leaving the total wave function in the usual 
gaussian-like solution. 

Note that in Figs, [l] and [2] the vortices appear in the 
"— " component. We have verified that, depending in the 
arbitrary initial conditions of the numerical method, the 
vortices may appear in any of them. As mentioned above, 
we believe there is an spontaneous symmetry-breaking 
mechanism acting here. However, as shown in Fig. |3| 
if the longitudinal magnetic field Bz — BiZq is different 
from zero, for instance zq = ±5, the vortices always ap- 
pear in the "±" component respectively. 

Fig. [4] is the 3D version of Fig. [l] We find that 
although the clouds have spherical symmetry, inherited 
from the trapping potential, the vortices are of cylindrical 
nature. This is because the nucleation points (where the 
magnetic field becomes zero) are actually filaments that 
cross the full BEC cloud. 

Fig. [5] shows single / = 2 vortices in a 2D BEC. To 
produce them, let / = 2 in Eqs.([2| and ([3|. Again, the 
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FIG. 1. Vortex of charge / = 1 at the origin fo = with 
J = 1/2. We show density plots, (a) component, (b) "-" 
component, and (c) total density; circulation fields, (d) 
component, and (e) "-" component; total circulation values 
as function of the radius of integrations circuit. Parameters 
are: g = 4000, k = 1.0 and zq = (i.e. no longitudinal field). 
The chemical potencial is /i = 29.9594 ± 0.0001, and the total 



relative spin populations are 



<Ar, > 



0.51 and 



<Ar_> 



: 0.49. 



The vortex is in " component and has circulation C/27r = 
-1. 



vortices appear at the chosen vanishing points of the mag- 
netic field. However, the vortices have charge / = 2, 
namely, the circulation gives C/27r = ±2. 

Results for J = 1 are reported in Fig. |6] The simplest 
case corresponds to a single / = 1 vortex at the origin. 
As predicted in Section II, the a = 1 and a = —1 compo- 
nents show single vortices with opposite sign circulations 
of charge 1; the a = component does not present a 
vortex. In this case, our method does not give rise to 
other charge distributions within the spin components; 
however, some may be found in the literature [25j . 
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FIG. 2. Vortex of charge / = 1 at the chosen place fo = (4, 2), 
with J — 1/2. We show density plots, (a) component, 
(b) "-" component, and (c) total density; circulation fields, 
(d) component, and (e) '-" component; total circulation 
values as function of the radius of integrations circuit. Pa- 
rameters are: g = 8000, At = 0.5 and = 0. The chemical 
potencial is /i = 46.5984 zh 0.0002, and the relative spin pop- 
ulation are 0.51 and 0.49. (b) The vortex 
is in " component and has circulation C/27r = —1. 



B. Multiple / = 1 vortices in 2D BEC 




FIG. 3. Vortices of charge / = 1 at the origin ro = 0, with J — 
1/2. (a,b) zq — 1, the vortex appears in the + component, 
with positive charge +1; chemical potencial is /i = 29.8472 ± 
0.0001 and the populations are — 0.37 and — 

0.63. (c,d) zq — —1, the location, population and circulation 
of the vortex is reversed; chemical potencial is /i = 29.8562 ± 
0.0006, populations are = 0.63 and = 0.37 




FIG. 4. Vortex of charge / = 1 at the origin fo = 0, with 
J —1/2 and zo — 0, for a 3D BEC. The vortex is in the " 
component with negative circulation. Chemical potencial is 
ji — 15.48 ±0.01; populations are — 0.53 and — 

0.47. 



It was reported in the previous subsection how external 
fields of the form given by Eqs.([2| and (|3| nucleate single 
vortices at the desired locations. We now turn to the 
possibility of nucleating arbitrary arrays of vortices. 

The first case is shown in Fig. [7| where we study a 
J = 1/2 2D BEC. The external magnetic field is a su- 
perposition of fields generated by eight "infinite" wires 
(with currents of different magnitudes and signs) placed 
outside the BEC cloud such that V • ^ = and V x 5 = 
inside the cloud region. A single wire located at fo has a 



field, 



B{x,y) = I 



{x - XQ)ey - {y- yo)e:, 
{x - xo)^ + - yo)'^ 



(14) 



where / is the current in dimensionless units, and is 
the unitary vector in direction ^. 

With the appropriate choice of the currents / (see the 
figure caption for details) , we were able to generate three 
spatial points where the field vanishes linearly, namely 
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FIG. 5. Vortex of charge / = 2 at the origin ro = 0, with J = 
1/2, ^0 = 0, ^ = 1000 and k = 0.2. The vortex is in the 
component with negative circulation C/27r = —2. Chemical 

<iV+> 



potential is /i = 13.9 ±0.3 
^=0.4. 



Populations are 



0.58 and 



/ = 1 "zeroes". This should generate three vortices at 
those locations and, as shown in Fig. [7| we do find 3 
vortices with / = ±1 charge. For the case = 0, two 
vortices appear in one component and one in the other. 
Nonetheless, in some cases the three vortices appear in 
only one component. Again, we have no explanation at 
this moment for these occurrences and we appeal once 
more to a symmetry breaking mechanism. However, if 
7^ 0, the three vortices appear (not shown) at the 
component corresponding to the sign of Zq. 

Fig. [8] shows a regular lattice of vortices. These were 
generated by a magnetic field of the form, 

B{x^ y) = B (sin(7rx/A) cosijiy / \)ex — sin(7r?//A) cos(7rx/A)6 

(15) 

The field vanishes linearly at the points {x^y) = 
(±nA, ±mA) with n and m both integers or both half 
integers; that is, those points are "zeroes" of dipolar or- 



FIG. 6. Vortices of charge / = 1 at the origin ro = 0, with 
J— I^zq— ^^g — 8000 and n — 0.5. Chemical potential is 
/i = 47.1020 zb .0003; populations are = 0.25, = 

0.25 and = 0.50 . The vortices appear at the mj — ±.1 

(a,b,d,e). There is no vortex at mj — (c). The vortex 
at mj = +1 has positive circulation C/2ti — +1, while the 
circulation is negative for mj — —1 (f) . 



der / — 1. 

Another interesting feature of this lattice is that vor- 
tices of charge / = ±1 are generated within the same 
spin component. In the zeroes of the magnetic field that 
correspond to both m and n integers, vortices of charge 
/ = +1 are nucleated in the "+" spin component and 
with charge / = — 1 in the "— " spin component. In the 
same fashion, in those zeroes corresponding to m and n 
both half integers, vortices of charge / = — 1 are nucleated 
in the "+" spin component and with charge / = +1 in 
the "— " spin component. Fig. [8] illustrates a vortex lat- 
^ ^ice of charge / = ±1 formed in the "+" spin component, 
following the rules just described. 

Such a behavior can be readily understood using the 
tools from Section II. Expanding the magnetic field 
around each of its zeroes, we note that the behavior of 
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FIG. 7. Three vortices of charge / = 1 at the chosen loca- 
tions ri — (3.6, —6.1), r2 — (—6.6, 0) and rs — (3.6, 6.1), with 
J = 1/2, = and n — 2.0. We show density plots, (a) 
component, (b) "-" component, and (c) total density; circula- 
tion fields, (d) component, and (e) '-" component; total 
circulation values as function of the radius of integrations cir- 
cuit. Note that the vortex in the "+" component has positive 
+ 1 circulation, while those in the " component each has 
-1 circulation. Chemical potencial is /i = 53.4084 ± 0.0002; 

<iV+> 



population are 



0.525 and 



: 0.475. 



the magnetic term in Eq.Q depends on the nature of the 
zeores. The magnetic term is 











(16) 



where is a constant. The first case corresponds to 
both n and m integers, while the second one to the half 
integer case. Again, the lattice may appear part in one 
spin component and part in another, unless a longitudi- 
nal magnetic field is applied. An "undesired" character- 
istic of this lattice is that V x 5 is not zero in the whole 
region occupied by the BEC cloud. That is, one needs 
"sources" in the region occupied by the condensate. Nev- 
ertheless, this case illustrates the possibility of generating 



FIG. 8. Vortex lattices of charge / = lforaJ=l/2 BEC, see 
text for magnetic field. In all cases g — 8000 n — (a,c) = 
showing the spontaneous symmetry breaking mechanism 
with chemical potencial \i — 49.5221 ±0.0003 and populations 
= 0.62 and = 0.38. (b,d) = +5 while (c,f) 

are for = —5 with chemical potencial \i — 40.3934 di 0.0006 
and \i — 40.3937 ± 0.0006 respectively. Populations in (b,d) 



are 
are 



N 
<iV+> 



0.004 and 
0.996 and 



N 
<N_> 



= 0.996; for (c,f) populations 
= 0.004 



vortices at the desired locations. We leave the tailoring 
of the corresponding fields to the ingenuity of the exper- 
imental researchers. This case can also be realized in a 
3D BEC. 



IV. FINAL REMARKS 

As we have seen, in multicomponent BEC clouds, the 
presence of transversal magnetic fields with vanishing 
values at given spatial locations, nucleate vortices at 
those very points. The vortices can thus be placed at 
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0.0 0.1 0.2 0.3 0.4 0.6 



xlO" 



FIG. 9. Vortex lattice corresponding to (b) in Fig. |8] The 
circulation field is shown in black arrows. One can observe 
the alternating character of the vortex topological charge. 



any location "on demand". We have verified this re- 
sult for many tailored magnetic fields. In addition, we 
have verified the further prediction that the total density 
p{x^y) = + \il)-{x^y)\^ does not show any ev- 

idence of the vortices and that its form is the "simple" 
solution of GP in a harmonic potential. 

It is interesting to contrast the present form of the 
appearance of the vortices with that of a rotating super- 
fluid. In the rotating frame, the (one-component) GP 
equation looks like [2 , 



2m 



V'^<^{r) + {n-L)<^{r)+g |$(r)r $(f) -/i$(r) = 0, 

(17) 

where Q = Qz is the angular velocity of rotation and L 
is the angular momentum operator. 



L 



:fx V. 



By a vector identity one can rewrite. 



(18) 



(19) 



and, therefore, one can see that the rotation term is 
equivalent to introducing a vector potential A = Ct x f 
corresponding to a longitudinal uniform "magnetic field" 
B = Q {in appropriate units). 

In its celebrated paper, Abrikosov[8 showed that if the 
above equations hold, a lattice of vortices, of the type 
described in Section I, will be created by the presence of 
such a magnetic or appropriate gauge field. It does ap- 
pear, therefore, that vortices are created by the coupling 
of the orbital angular momentum L with the external 
field ft, i.e. ^1 • L. The fact that the vortices appear in a 
lattice is "incidental" due to the global structure of the 
GP equation in the presence of a uniform longitudinal 
magnetic field. 



In the case described in this work, vortices also appear 
because of a magnetic angular momentum coupled to an 
external gauge field. It does appear, however, that in 
the present case it is the intrinsic or spin angular mo- 
mentum the main player instead of the orbital angular 
momentum. 

We believe that the above two cases may not be as 
different as they may appear. This assertion may be 
based on the observation that GP equation near a zero 
of the magnetic field can be transformed as one with 
an artificial gauge field coupled to the orbital angular 
momentum. Let us see this briefly. For simplicity, take 
the case J=l/2, zo = and g^fd = g^ajs- The coupled 
equations ([2| may be written as 



2m 
where , 



^ = 0. 

(20) 

(21) 



and = + By applying the unitary trans- 

formation 



U = g-*0crz/2g-i7rcry/4gi(/)cr^/2 



(22) 



with tan(/) = y/x, equation (20) now becomes [ 

" 1 

2m ' 

^Ve:,t{r) ^ 9\^\^ ^ i^ra, - ^ ^ = 0, (23) 



-{-ihV - AxGx - AyGy - AzCFzY 



with ^ the transformed state. The induced artificial 
gauge potentials may be easily calculated [33]. In par- 
ticular, the diagonal potential Az is. 



A.= 



n 

Yr' 



(24) 



which corresponds to a nonuniform singular magnetic 
field, actually, to that of a magnetic monopole. Thus, 
the vortex nucleation in the rotating condensate ap- 
pears to be equivalent to ours via an artificial mag- 
netic monopole field. This, in turn, is intimately related 
to Berry phases |45j. Therefore, the present scheme 
may also be seen from the perspective of artificial gauge 
fields [T6l [T7| [32 1 146 } [5T] . although it should be clear that 
we do not assume adiabaticity of the diagonal terms in 
([23|. 

We believe the present results can be simply realized in 
actual experimental conditions. In particular, for a ^^Rb 
BEG there are two cases that can be simply matched. 
One is the F = 1 case that naturally corresponds to the 
J = 1 situation here discussed. The other, the case F = 2 
may be tailored to have only the components mi = 1 
and m2 = 2 corresponding to J = 1/2. The external 
magnetic fields may be produced by single- wire magnetic 



9 



fields, as given by (14). However, an additional non- 



uniform small magnetic field in the z direction should 
be used to separate the spin components such that the 
phenomena may be observed. 

The present scheme may also be used as a route to 
quantum turbulence. Recently, Bagnato et al.[39 found 
a way to stir a ^^Rb condensate by means of time de- 
pendent magnetic fields, achieving a state of quantum 
turbulence. Although we are not asserting that our pro- 



cedure explains those experimental results, it does appear 
that once vortices are nucleated by transversal magnetic 
fields, the vortices may be stirred by simply changing in 
time the positions of the zeroes of the magnetic field. It 
seems evident that a strong stirring of the vortices may 
lead to a quantum turbulent state. The study is under 
way and will be reported elsewhere. 

We thank support from grant DGAPA UNAM 
IN108812. R.Z.Z. and M.L.H acknowledge support from 
CONACYT, Mexico. 
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